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It has been recently shown that a kinetic theory can be modified to incorporate 
triangle anomalies and the chiral magnetic effect by taking into account the Berry 
curvature flux through the Fermi surface. We show how such a kinetic theory can 
be derived from underlying quantum field theories. Using the new kinetic theory, 
we also compute the parity-odd correlation function that is found to be identical to 
the result in the perturbation theory in the next-to-leading-order hard dense loop 
approximation. 
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I. INTRODUCTION 

Kinetic theory [1] has wide applications in condensed matter physics, nuclear physics, 
astrophysics and cosmology. There is, however, an important deficiency of the conventional 
relativistic kinetic framework: it misses the effect of triangle anomalies [2, 3] — an important 
feature of relativistic quantum field theories. Recently it has been shown in Ref. [4] 1 that 
a kinetic theory for Fermi liquids can be modified to include such anomalous effects by 
taking into account the Berry phase and Berry curvature [12] - the notions extensively 
studied and widely applied in condensed matter physics [13]. It was shown that not only 
the form of the transport equation but also the definition of the particle number current 
must be modified when the Berry curvature has a nonzero flux through the Fermi surface. A 
consequence of this modification is generation of parity-violating and dissipationless current 
in the presence of magnetic field, called the chiral magnetic effect [14-17]. This had been 
previously found in the perturbation theory [14, 17] and the gauge/gravity duality [18, 19] 
and was incorporated in the framework of hydrodynamics [20] (see also Refs. [21, 22] for more 
recent developments). The chiral magnetic effect may have been experimentally observed 
in relativistic heavy ion collisions [17, 23] and is potentially observable in Weyl semimetals 
which possess band-touching points [24-26]. 

On the other hand, one should be able to derive the kinetic theories from the underlying 
quantum field theories by following the standard procedure: starting from the equations of 
motion for the two-point function {ip(x)ip^ (y)) and performing a derivative expansion for its 
gauge-covariant Wigner transform, one arrives at the Vlasov equation (see, e.g., Ref. [27] 
for a review). So far, Berry curvature corrections to the relativistic kinetic theory have been 
ignored in the field theoretic derivation. Also, microscopic origin of the modification to the 
particle number current is not yet clear. 

In this paper, we microscopically derive the kinetic theory with Berry curvature correc- 
tions from underlying quantum field theories. 2 For concreteness, we consider the system of 
relativistic chiral fermions at finite chemical potential /i (which is known to have a nonzero 
Berry curvature flux). Our starting point is the high density effective theory [30, 31] that 
describes the physics near the Fermi surface of chiral fermions. In this effective theory, one 
decomposes two-component chiral fermions into single-component particles ip + with positive 
energy E — \p\ — fj, and antiparticles ip_ with negative energy E = —\p\ — fi. Then one 
usually concentrates on the former with E ~ for |p| ~ fi, while neglecting the latter with 
E ~ — 2/x; picking up only ifj + degrees of freedom leads to the conventional Vlasov equation. 
As we shall demonstrate in this paper, however, if one carefully integrates out ip- degrees of 
freedom, Berry curvature corrections emerge in the kinetic theory from the mixing between 
ip + and ip- (or and i])-). Modification to Liouville's theorem on the phase space known 
in the condensed matter literature [32, 33] and modification to the current found in Ref. [4] 
can be naturally understood from this deliberate integrating out procedure [see Eqs. (72) 

1 See Refs. [5-9] for further investigations and applications. See also Refs. [10, 11] for different approaches 

to derive kinetic equations with triangle anomalies without referring to the Berry curvature. 

2 See Refs. [28, 29] for related attempts. 
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and (77)]. Apparently, the essential ingredient in this field theoretic argument to lead to 
Berry curvature corrections is a Fermi surface of chiral fermions. 

We also compute the parity-violating correlation function, using the kinetic theory with 
Berry curvature corrections. In the case of the conventional Vlasov equation, it is known 
that the parity-even correlation function computed in the kinetic theory coincides with the 
one in the perturbation theory under the hard dense loop approximation [27, 34]. In this 
paper, we will see that the parity-odd correlation function derived from the new kinetic 
theory is equivalent to the result in the perturbation theory beyond the leading-order hard 
dense loop approximation. 

The paper is organized as follows. In Sec. II, we review the kinetic theory with Berry 
curvature corrections. We also derive a new relation for the spin magnetic moment of 
quasiparticles in Fermi liquids. In Sec. Ill, we derive the new kinetic theory starting from 
quantum field theories. In Sec. IV, we compute the parity-violating correlation functions 
using both the new kinetic theory and perturbation theory, and confirm their agreement. 
Section V is devoted to our conclusions. 

Throughout the paper, we consider sufficiently low temperature regime T <^ fi where the 
Fermi surface is well-defined. We also concentrate on the collisionless limit of the kinetic 
theory. 



In this section, we review the kinetic theory in the presence of the Berry curvature which 
exhibits triangle anomalies [4] (see also Refs. [7, 9]) and provide the proper definitions of 
particle number density and current. We also derive the dispersion relation of quasiparticles 
according to the constraints of Lorentz invariance. 



We first consider a single chiral fermion expressed by the two-component spinor u p sat- 
isfying the Weyl equation 



where the signs + and — correspond to right-handed and left-handed fermions, respectively. 
The two-component spinor described above has a nonzero Berry connection defined by [12] 



II. KINETIC THEORY WITH BERRY CURVATURE 



A. Berry curvature and Poisson brackets 




(1) 



p- 



(2) 



and a nonzero Berry curvature, 




(3) 
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where p = p/|p| is a unit vector. Equations (2) and (3) can be regarded as the fictitious 
vector potential and magnetic field in the momentum space. This fictitious magnetic field 
can be associated with the one from a "magnetic monopole" with the charge ±1/2 put in 
the center of the momentum space. As a result, the motion of chiral fermions is affected by 
the Berry curvature in the momentum space, in addition to the usual electromagnetic fields 
in the coordinate space. In particular, the effects of the Berry curvature work oppositely 
between right-handed and left-handed chiral fermions. 

Let us now consider the action of a single quasiparticle in the presence of the electromag- 
netic fields and Berry curvature [32, 33], 

S = J dt^x 1 + A^x)? - A\p)p i - e p (x) - A°(x)}. (4) 

Note that the quasiparticle energy e p is a function of x in general; indeed chiral fermions 
have the magnetic moment at finite chemical potential \x and their energy depends on the 
magnetic field B(x) [see Eq. (43) below]. The action (4) can be summarized in the following 
form, by combining space x and momentum p into a set of variables £ a (a — 1, • • • , 6), 

S = J dt[-u a (t)?-H(t)], (5) 

where H(£) = e p + A is the Hamiltonian. 

The equations of motion of the action (5) read 

oo ab i h = -d a H, (6) 

where u a b = d a 0Jb — dbU a and d a = d/d^ a . This equation can be rewritten as 

t = -u ab d b H, (7) 

where u ab = {u~ l ) ah is the inverse matrix of uj a b- Here we assume the existence of the inverse 
matrix, i.e., u = detu a b ^ 0. Equation (7) can be interpreted as 



dH 

once we define the Poisson brackets as 



{UZ b } = to ab . (9) 
The explicit forms of the Poisson brackets for the action (5) read [33] 

r„ „ i _ e ijkBk r -j _ ejjkQk r -j _ Qjj + ^jBj , , 

where B l = e l ^ k dA k /dx 3 ' . These Poisson brackets should be compared with the usual ones 
in the absence of the Berry curvature and electromagnetic fields: 

{Pi,Pj} = ®, {xi,Xj} = 0, {pi,x j } = 5 ij , (11) 

whose invariant phase space is <ip<ix/(27r) 3 . As a consequence of the modifications to the 
Poisson brackets above, the invariant phase space is modified to [32] 

dr = v^d£ = (l + B.fi p )|^, (12) 

where oj = det u a b- 
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B. Kinetic theory with Berry curvature and triangle anomalies 



Let us construct the collisionless kinetic theory incorporating the effects of the Berry 
curvature. If collisions between particles are negligible, each particle constitutes a closed 
subsystem. According to the Liouville's theorem which states that a volume element in the 
phase space does not change during its time evolution, the one-particle distribution function 
n(£) would obey dn/dt = 0. However, the invariant phase space is modified as Eq. (12) 
due to the Berry curvature, and the probability of finding a particle in the phase space is 
\/un(C,)dC,. As a result, we instead use the modified distribution function p(£) = v /cJn(£) 
that obeys the equation dp/dt = 0, or equivalently, 



P + d a (i a p) = o. 



Using Eq. (7), this reduces to 



h p - uj ab d b Hd a n p = 0. 



(13) 



(14) 



Setting H = e p + A , we can explicitly write down the kinetic equation (see also Refs. [7, 9, 
33]): 



Hp + 1 + B • tt r 



(E + v x B + (E-B)ri p ) 



dn 



^+(v + Extt p + (v.ft p )B) 



drip 
<9x 



= 0, 
(15) 



where v = de p /dp and E = E — <9e p /<9x. This is a general low-energy effective theory in the 
presence of Berry curvature corrections that describes the evolution of n p . Note that v is 
different from the unit vector p when the quasiparticle energy e p has the contribution from 
the magnetic moment [see Eq. (43) below] . If we turn off the Berry curvature and ignore the 
x-dependence of e p (i.e., £1 = and de p /dx. = 0), this reduces to the usual Vlasov equation. 
We now define the particle number density 

d 3 p 



n 



(2nf 



;i + B-n p )n p , 



(16) 



where the invariant phase space is modified according to Eq. (12). Multiplying the kinetic 
equation (15) by y/u, performing the integral over momentum p, and using Maxwell equa- 
tions, V • B = and 9 t B + V x E = 0, we obtain the following identity [4]: 



d t n + V • j = 



where 



J 



d 3 p 
(2tt) 3 



drip 
' dp 



+ n 



d 3 p 
(27p 



drip 
dp 



drip 
dp 



E B, 



e p B + e p ri p x — 

ox 



+ E x cr, 



is identified with the current and cr is defined as 

d 3 p 



f d 6 i 

a = JW; 



3^ p n p . 



(17) 



(18) 



(19) 
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In Eq. (17), we observe that the particle number of chiral fermions is no longer conserved 
when we turn on both electric and magnetic fields. This can be evaluated, by integration by 
part and using Yp • fl p = around the Fermi surface, n p — 1 deep inside the Fermi surface 
and n p = far outside the Fermi surface, as 



9 t n + V.j = ±^E.B, 



(20) 



for right-handed and left-handed fermions, respectively. This is exactly the equation of 
triangle anomalies in relativistic quantum field theories [2, 3], which holds independently of 
interactions. 

The first term in Eq. (18) is the usual particle number current of the kinetic theory, while 
the remaining terms are the Berry curvature corrections. The same form of the current can 
be obtained in the Hamiltonian formalism using the commutation relations postulated in 
Ref. [4]. The final term in Eq. (18) is the anomalous Hall current, which vanishes for a 
spherically symmetric distribution function at rest. In this case, the current is 



f*2 

J (27T 



d 3 p 

Y 3 



+ n 



On 

e p B + e p Q p x 



(21) 



~ p dp ' V"" P dp 

At this moment, microscopic origins of the Berry curvature corrections to the particle number 
density and current in Eqs. (16) and (21) are not so clear. In Sec. Ill, microscopic meanings 
of these corrections will be clarified in the field theoretic language. 

It should be remarked that there is an ambiguity to define the number current from 
the continuity equation because j = j + V x a with any vector a is also a solution to the 
continuity equation. In order to fix this ambiguity, we look at the energy and momentum 
conservations. Define the energy density and the momentum density, 

d 3 p . . if d 3 p 



I 



1 + B • ft p )e p n p , ir l 



I 



;i + b • n p yn p . 



(22) 



(2tt) 3V ' "P^P'-P' J ( 27r )3' 

Multiplying Eq. (15) by e p ^Juj and p l y/u, and performing the integral over momentum p, 
we have 



d t e + 



d 3 p 



e p F = 0, 



d 3 p 



p l F = 0, 



(23) 



(2tt) 3 p " i " ' J (2tt) 3J 

where F is the piece in the square brackets of Eq. (15). The above equations can be 
respectively interpreted as the energy and momentum conservation laws 



E i j i , dpT* = nE i + e ijk j j B k , 



where 



rpOi _ 



/ 
/ 



d 3 p 
(2^)3 
d 3 p 

W) 



(S ij + Fff) 



e n dn T 



IP 



2 dpi 



+ e 



ijk Ppj p 



dx k 



dp k 



E l n p + e p 



drip 
dx l 



(24) 
(25) 

5 ij e, (26) 
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which indicates that j l is the genuine current. 

Alternatively, this ambiguity is avoided if we use the definition of the current [7, 9] 




By using Eq. (7), one can actually check that this current is equal to Eq. (18). 

In equilibrium where n p is homogeneous, the first and third terms in the right-hand side 
of Eq. (21) vanish, while the second term is nonvanishing. Using e p = /x at the Fermi surface 
with /i the chemical potential, we find 

This is the relation of the chiral magnetic effect [14, 15, 17]: an equilibrium current induced 
in the direction of the magnetic field for chiral fermions at finite chemical potential \x. 



C. Lorentz invariance in Fermi liquids 



Here we 
dau's Fermi 



consider the consequences of Lorentz invariance in a system described by Lan- 
i liquid theory [35]. The constraint due to Lorentz invariance is that the energy 
flux is equal to the momentum density, T° l = ir\ From Eqs. (22) and (25), this condition 
in the homogeneous system becomes 

" / Wf {Sii + **4w = I W? (1 + B ' n * )p ' n - (29) 

We vary both hand sides of Eq. (29) as n p = n p + 5n p and e p = e p + 5e p , where 

Se p = J -J^ (1 + B • n q )/(p, q)5n q , (30) 

with /(p, q) being some function characterizing the interactions among quasiparticles, called 
the Landau parameters. By performing the integral by parts, we have 

/ + *H) gigr** " 4*pg] =Jw? il + B - n " )p ' Sn "- (31) 

Using Eq. (30) and renaming variables q -B- p, the second part of the left-hand side reduces 
to 

-/(& < x + B - "J + B ' n ' )/(p ' Sn - (32) 

Because Sn p is arbitrary, we have the relation: 

<9e° f H 3 n dn° 

+ B HK-^ - (i + b • fi p ) J ^L(s- + b^Kp, q ) £q ^ = (i + b ■ n p y. 

(33) 
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To proceed, we take an ansatz 



4 = v f(P ~ Pf) + ^ + t(p) b • P, 



(34) 



to the linear order in B, with some constants Vf, pf, fi, and 7 being a scalar function 
of p = |p|. Note that the Fermi velocity defined by de^/dp is B-dependent. Note also 
that Landau parameters are functions of B; from the property /(p, q) = /(q, p), Landau 
parameters are composed of two parts (to the linear order in B): 



/(p,q) = /Vq)+ B 'f 2 + ' ) /Vq), 

zp f 



(35) 



where f A ' B (p, q) are independent of B and can be expanded by the Legendre functions as 

00 

f A ' B (p,ci) = J2ff' Bp i(^e), (36) 



1=0 



where 9 is the angle between p and q both taken on the Fermi surface. 

We now evaluate both hand sides of Eq. (33) to the linear order in B. Substituting 
Eqs. (34) and (35), replacing p by pf, and performing the angular integral (note also that 
rip has the B- dependence), we have 



liv f [l + -Ff ) f 4 



Pf (v f i + H - ^v f F[ A 



+ 



2p) 



3 1 3 1 5 2 



2p 2 f 



1 + F q a + \f q b -^F 2 b 



+ 7 



B l 



(B • p)p l 
(37) 



where 7 = 7(p/), 7' = §^l{Pf)i an d we defined 



/ 



dq i 



(B q) = -f[ (B ■ p)p\ 



(38) 



In order to satisfy Eq. (37) for any B and p, we must have 



7 



2 P j 



Pf (v f l + H - \^v f F[ A ^ + || 



+ 0+ 3° 15 2 



(39) 
(40) 
(41) 



While Eq. (39) is the relation obtained by Baym and Chin [36], Eq. (40) is a new relation 
for the anomalous spin magnetic moment. A constraint from the gauge invariance on the 
anomalous angular magnetic moment in Fermi liquids was originally given by Migdal [37] 
and was studied in details in Ref. [38]. Here we have provided the new constraint on the 
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anomalous spin magnetic moment from the viewpoint of the Berry curvature together with 
the Lorentz invariance. 

In particular, in the noninteracting limit where /(p, q) are turned off and pf = /i, we 
obtain 

^ = 1 > TM = -2^, 7 ' (/i) = ^ (42) 
A solution to satisfy these relations is taken as 7(2?) = — l/(2p 2 ). In this case, we have 

We shall see in Sec. Ill [Eqs. (52) and (66)], based on the microscopic quantum field theories, 
that this is actually the dispersion relation of chiral fermions near the Fermi surface. 



III. FROM QUANTUM FIELD THEORIES TO KINETIC THEORY WITH 

BERRY CURVATURE 

In this section, we derive the kinetic theory constructed in Sec. II B from the microscopic 
quantum field theories: the kinetic equation (15) and the modified number density (16) and 
current (18) are reproduced microscopically. The resultant kinetic theory exhibits triangle 
anomalies and the chiral magnetic effect. 

Our procedure is as follows: we first consider the high density effective theory [30, 31], 
which is an effective field theory valid near the Fermi surface. The expansion parameter of 
the theory is l/fj, where / is the residual momentum measured from the Fermi surface. We 
then derive the kinetic theory, by performing the derivative expansion for the equations of 
motion of the Wigner function defined in the high density effective theory. The expansion 
parameter here is the slowly varying disturbances dx taken to be much smaller than the 
chemical potential /j, and the gauge field A^. Note that our procedure does not rely on 
the expansion in terms of the coupling constant, and hence, it is applicable even when the 
interactions are strong as long as the notion of quasiparticles is well-defined. 

In this and following sections, we consider the theory with right-handed chiral fermions. 



A. High density effective theory 

We first review the derivation of the high density effective theory [30, 31]. We start with 
the Lagrangian for right-handed fermions, 

£ = ^(ifl + fi)ij, (44) 

where Ip = a^D^ with = <9 M + iA^ and <j fl = (1, <x), and A^ is the external background 
field. In order to focus on the particles near the Fermi surface, we decompose the energy 
and momentum of a particle (or a hole) near the Fermi surface as p° = /x + 1° and p = /iv + 1 
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with 1°, |1| <C /i, where v is a unit vector which specifies a direction to a point on the Fermi 
surface. The momentum can be shifted by /xv by performing a Fourier transformation 

^(x) = J2^ V ' X Mx), (45) 

V 

where summation is taken over v. The matrix a ■ v in the momentum space can be diago- 
nalized by using the projectors as 

= P±(v)1>v, P±(v) = l± l' W - (46) 

In the absence of the external electromagnetic field, ip± satisfies the eigenvalue equations, 
(<r • v)V>± = ±^±. 

In terms of ip± v , Eq. (44) reduces to 

V> f (^ + /i)^ = Yty+v iv ' + ^-«( 2 A* + ™ ■ D )^-v + {^\ v %V>i^-v + h.c.)], (47) 

where = (l,v), ^ = (1,-v), < = (0, a - v(v • o-)), D\ = (0, D - v(v • D)), and 
fl ± = cr^D^ = a^D^. From Eq. (47) we obtain the equation of motion for ip- v , 

(2/i + iv ■ D)ij_ v + ifl ± ip +v = 0. (48) 

By integrating out ip- v using Eq. (48), the effective Lagrangian in terms of ip + can be written 
down order by order in l//x: 

£eft = E £(n) ' £(n) = E (49) 

n v 

where denotes the effective Lagrangian of the n-th order in lf/i (n — 0, 1, 2, • • • ). The 
explicit expressions for (n = 0, 1, 2) are 

2><°> = it, • A ^ (1) = |i, © (2) = "^JM* • Wl- (50) 

Using = D\ + B • cr and xjj\. v crijj +v — w+v^+v, the dispersion relation near the Fermi 
surface to 0(l/fi) reads 

Taking into account |p| = /x + + l\/{2ji) + 0(l//x 2 ), it can be rewritten as 

B v ^ / 1 

where p = |p|. This indeed agrees with Eq. (43) up to this order. (The agreement will be 
shown to higher order in 1/// in the next subsection.) The second term in Eq. (52) originates 
from the magnetic moment of chiral fermions at finite fj,. This is similar in structure to the 
Pauli equation that describes the magnetic moment of massive Dirac fermions in the vacuum; 
the Pauli equation can be obtained by expanding the massive Dirac equation in 1/m, where 
m is the mass of Dirac fermions. 



Zv=P- — + 0[-), (52) 
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B. Kinetic theory via derivative expansion 

We construct the kinetic theory based on the effective theory (49) on a patch indicated 
by a unit vector v. We consider the Dirac operator to the second order in 1/fi, T> = 
p(o) _|_ _|_ -p( 2 ) ) anc i introduce a two-point function for ip v , 

G v (x,y) = (,f> v (x)if>l(y)). (53) 
The function G v (x,y) satisfies equations of motion together with projection conditions, 

V x G v (x, y) = 0, G v (x, y)V\ = 0, (54) 
P-(y)G v (x,y) = 0, G v (x,y)P-(y) = 0. (55) 

In thermal equilibrium where the system is homogeneous, G v depends only on the relative 
coordinate s = x — y. We are interested in the small deviation from the equilibrium where 
G v depends both on x and y. It is thus useful to change the coordinates from (x, y) to the 
center-of-mass and relative coordinates (X, s) defined by 

and consider the derivative expansion with respect to X. In order to derive a quantum 
analogue of the classical distribution function, we perform the Wigner transformation: 

G V (X, l) = J d A s e ihs G v (X + |, X - , (57) 

where Z M is the residual four-momentum. Unlike G v (x,y), however, G v (X,l) is not gauge 
covariant. We will use the gauge covariant definition instead: 

G V (X, l) = J d 4 s e u - s U (X,X+ S -)G V (X+ S -,X- S -)U(X- S -,X), (58) 



where 



U (x, y) = Pexp 



dx^A^x) 



(59) 



is the Wilson line. The symbol P the path-ordering along the path 7 from x to y. For 
simplicity, G v (X,l) is renamed G v (X,l) in what follows. 

In constructing the kinetic theory, we consider the slowly varying disturbances and per- 
form a gradient expansion in terms of dx- To this end, we assume the following counting 
scheme: dx = 0(ei), d s = 0(e 2 ), = 0(e 3 ), and F^ u = 0(eie 3 ). Here dx and d s are 
the derivatives with respect to X and s, and q (i = 1,2,3) are independent expansion pa- 
rameters which satisfy the conditions, e\ 62,3 1. The condition 1 {% — 1,2,3) is 
necessary for the derivative expansion in the high density effective theory while e\ ^ e 2j 3 is 
necessary for the derivative expansion in the kinetic theory. In order to take into account 
triangle anomalies, we consider the kinetic theory to 0{e\e^) [or the occupation number n p 
to 0(eie§)]. 
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For simplicity, in this subsection we consider the homogeneous system where d x F^ u = 
and dfrip = (n p is the distribution function which will be defined below). This is sufficient 
for our purpose to understand the microscopic origin of the Berry curvature corrections. The 
generalization to the inhomogeneous case should be straightforward. 

Consider the equations, T> x G v (x,y) ± G v (x,y)V^ y = 0. We expand them in terms of dx 
to the second order and perform the Wigner transformation. The Wigner transform of the 
equations can be written down order by order, 

I { ± ] = J ^- 4 e lls (V^G V ± G B 2>Wt), (60 ) 
for n — 0, 1, 2. The expansion of the gauge field in dx reads, 



A,(x) » A,(X) + \{s ■ 8 X )A,{X) + l{s ■ d x fA,{X). 



(61) 



Combined with the contributions from the Wilson loop in Eq. (58), all the terms involving 
the gauge field are expressed by the gauge invariant field strength F^ v at the end. Then the 
third term in the right-hand side of Eq. (61) will be irrelevant eventually when dfF^ = 0. 
Renaming the kinetic residual momentum 1^ = 1^ — A^ as Z M , we have 



/| 0) = 2(/ -/||)G„ 



r(0) 



= iv»{g lia d t -F llv dr l )G v , 
(i) ' 



—l±(dfiodt — F^ u di)G v , 



= -(-ll + B-v)G v , 
1+ = ^Wl - B • v) + B • 1 ± + (E x 1) • v]G v , 



r(2) 



2/i 2 



1 



2 ¥1 - i^l - B • V K - e ijk v k v a F ia gM 



(g^odt — F^ u di)G v , 



(62a) 
(62b) 

(62c) 
(62d) 



where d t is the derivative with respect to the residual momentum /, and 1° = l\\ +0(/ 2 //i) is 
used in Eqs. (62c) and (62d). 

From the equation V x G v (x,y) + G v (x,y)V y = 0, we obtain the on-shell condition. Con- 
sidering the projection conditions (55), G v can be written as 



G v = 2irP + (v)6 No- h - 



Zi-B-v /||(/i-B-v) + B-lj 



+ 



2, ■ v < 63 > 

where ni(X) is the distribution function expressed by the residual momentum /. Recalling 



p — n + l\ + 



I 2 



Z=? Zi. 



+ 



Bp 

2p 



B v B 1 

1 ; 

2/i 



2/x 2/i 2 

ZiiB-v 



2/i 2 



1 

+ 



(64) 
(65) 



the dispersion relation in the delta function of Eq. (63) is equivalent to the condition p° = e p 
with 



e P =P 



Bp 

2p ' 



(66) 
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which indeed coincides with Eq. (43). Accordingly, the distribution function rii can be 
replaced by n p in terms of the original momentum p. 

On the other hand, from the equation T> x G v (x,y) — G v (x,y)V^ y = 0, we obtain the 
transport equation for ni(X). Using 



v + 

H 



2/i 2 



+ 0| .S 



the transport equation can also be expressed by the original momentum p. We end up with 
the transport equation: 



1 + 



B 



2/i 2 



:E + pxB) + (E-B) 



2/i 2 



dn p 
dp 



= 0. 



(68) 



Equation (68) indeed agrees with the homogeneous limit of the kinetic theory (15) to 0(efe§) 
if we identify fl p = p/(2p 2 ) at p — fj,. Therefore, we have found that the Berry curvature 
corrections in the kinetic theory emerge as the higher-order corrections in l//i to the usual 
Vlasov equation. Once the kinetic equation (68) is obtained, the relation of triangle anoma- 
lies (20) follows, as we have seen in Sec. II B. 



C. Particle number density and current 

Here we consider the particle number density and current for right-handed fermions with- 
out reference to the kinetic theory derived above. To see the chiral magnetic effect, we need 
to consider the number current to 0(6163) in the high density effective theory. Our discussion 
in this subsection is applicable to inhomogeneous electromagnetic fields. 

By definition the number density of right-handed fermions consists of four parts: 

= n ++ + n + _ + n_ + + n , (69) 

where n ++ is given by 

f d A p f d 3 p 

and n + _ = n_ + = because of the property of projectors, P + (v)P_(v) = 0. In order to 

express n in terms of ip+ v , we use Eq. (48) which relates ip- v to ip +v . Then n is given 

by 

Using the momentum p, we have in total 
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This is the number density including the Berry curvature correction, Eq. (16). 
Similarly the number current of right-handed fermions is decomposed as 

= j++ + j+- + j-+ + j— , 



(73) 



where j ++ is given by 



J++ 



d A p 



d 3 p 



(2tt) 4 v " j J (2tt) 3 
Using Eq. (48), summation of j + _ and j + can be written as 

d A p 



vn t . 



(74) 



2/dJ (2tt) 4 



tr[(£>i ± + i^ k v k Di)G v - G V {D^ - ie^v k D\% (75) 



while j ~ I//-* 2 is higher order in 1/fi and is negligible to the order under consideration. 

One can then rewrite Eq. (75) by changing the coordinates from (x, y) to the center-of- 
mass and relative coordinates (X, s) and performing the Wigner transformation in a gauge 
covariant way. One finds 



f + -+f- 



d 3 p 1 



(2vr) 3 2/x 

Putting them together and writing in terms of the momentum p, we arrive at 



(76) 



JR 



J (27T 



d 3 P _ 

) 3 



— B 

dp* p \2n dp 



p dn A c ijk ^ drip 



2/i dX k 



(77) 



where we used the dispersion relation (52). This is the same form as the current (21), includ- 
ing the chiral magnetic effect and the inhomogeneous term; the Berry curvature corrections 
in Eq. (21) microscopically originate from the mixing between ip + and ■?/>_. 



IV. CORRELATION FUNCTIONS 

In this section we compute the one-loop polarization tensor in the presence of chiral 
fermions at finite chemical potential // at zero temperature, using the perturbation theory 
and the kinetic theory constructed in Sec. II B. We confirm that both calculations give 
the same result, not only to the leading order but also to the next-to-leading order in \j \i. 
In particular, we find the kinetic theory with Berry curvature corrections reproduces the 
parity-odd polarization tensor beyond the leading-order hard dense loop approximation in 
the perturbation theory. 
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A. Perturbation theory 

The time-ordered Dirac fermion propagator at finite chemical potential /j, and zero tem- 
perature is given by 

S(x,y) = (Tj,{x)$(y)) = J - y )^^ + fa**^) 

-9(y° - x°)((3 p e^ x - y) + a p e ip(x - y) )] , (78) 

where a p = 9(p — /i), (3 P = 9(/i — po), a p = 1, and j3 p = 0. The one-loop polarization tensor 
in the presence of chiral fermions is then 

U^x-y) = ±tr[(l + >y 5 )>fS(x,y)>fS(y,x)]. (79) 
In the momentum space, it reads (see Ref. [39] for the case of Dirac fermions), 



d A x 



(2vr) 4 



e ikx W y (x) 



1 r _d*i 

2 J (2tt 



d 3 q_J_J_ 



T^(p,q) 
+T^(p,q) 
+T^(p,q) 
+T^(p,q) 



Po - qo - k - ie Po~qo-k + irj 

OL p OLq PpPq 



Po + qo - h - ie Po + qo~ h + if] 

Po + qo + ko-irj p + q + k + ie 

CX p j3q (3 p C\q 

Po - qo + k - it] Po - qo + h + ie 



,(80) 



where p = (p ,p), p = (po, -p), q = (<7o,q), q = (qo, -p), Po = |p|, go = |q|, p = q + k, and 
T^^p, q) = tr[(l + 75)7^7^ An infinitesimal quantity 77 takes rj = e for the time-ordered 
function IT^, and rj = — e for the retarded function n^. Substituting explicit expressions 
of the distribution functions, we find the /z-dependent part of the retarded function, 



n £ (t) 1 t*< 1 1 



2 J (2vr) 3 2po2g, 







[Ofr - a) - Obi - po)]- l ' : {P ' ! ' ] 



Po- qo~ h- ie 



Po + go ~ ^0 - «e Po + go + ^0 + «e 

where il-t and denote parity-even and odd parts, 



(81) 



Tnp,g)=tr( 7 Wg% (82) 
T^(p,g)=tr( 757 ^ 7 ^. (83) 

The parity-even part, n+", is the leading contribution in the hard dense loop approximation, 
while the parity-odd part, IT^, is suppressed compared with n^" by a factor of kj\x. 
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For completeness, let us first recall the computation of n+". Under the hard dense loop 
approximation (where k , |k| <C //), we end up with [34] 



= -\J 



''''-Mf, q ) (r'V • /•''/•" 2co 



2ir 2 



(2k) 



S^°S u0 - co 



v ■ k + ie 



4:71 v ■ k + ie 



(84) 



where v = (1, v), v = (1, — v) with v = q/|q|. 

Let us now turn to the parity-odd part Using tr(7 5 7 M 7 a 7 iy 7 /3 ) = — 4ie mul3 , we 

have 



nr(fc) = -- 



2 J (2tt) 3 4|q| 



A- fj.ua/3u Vo \ n \ 

k ' ^ ~ |q|) - T - ^ - |q|)2i^(fc a ^ + k a v p ) 



k ■ v + ie 



(85) 



The second term is vanishing while the first term remains nonzero only when (a, (3) = (0, k) 
and (k,0). Collecting both contributions, we obtain the expression: 

dv <jjv k — k k 



iri(Jfe) 



JL 

Att 2 



ie ijk k k + ie ijk 



to 



Att k ■ v + ie 



(86) 



where k denote the spatial indices [IT_^(A;) is vanishing otherwise]. Performing the an- 
gular integration, we finally arrive at 

,2 



m(k) 



4tt 2 



ie ijk k k 1 



CO 



[l-coL(k)], 



where 



L(k) 



In 



CO + 


|k| 


CO — 


|k| 



(87) 



(88) 



Equation (87) reduces to a simple form in the static or long wavelength limit: 

; ie ijk k k [to < |k|) 



iTi(k) 



_fU 

4^2 



' ie ijk k k (to > Ikl 



(89) 



127T 2 



That the latter, limiki/u,-^ lVi(k), is smaller than the former, linwiki-^o ^Vi(k), by a factor 
of 3 is consistent with the result of the "chiral magnetic conductivity" in Ref. [41]. 



B. Kinetic theory with Berry curvature 

We now compute the same retarded correlation function from the kinetic theory (15) 
constructed in Sec. II B, through the linear response theory: 



f(x) = J d% W\x-y)A v ( y \ 



(90) 
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or in the momentum space, 

f(k) = W v (k)A v (k). (91) 

Here we are interested in the current induced by a linear-order deviation of the gauge field 
Ap. For definiteness, we set up the following power counting scheme: A^ = 0(e) and 
d x = 0(5), where e and 5 are small and independent expansion parameters. Under this 
counting scheme, we compute the deviation of the distribution function Sn p and the current 
f to 0(e5). 

Remembering that dn p /dx. or de p /dx. can be nonvanishing at least when dn p /dx. = 0(e) 
or de p /dx. = 0(e) in Eq. (15), it is sufficient to consider the following kinetic equation of 
order 0(e5 2 ) [or n p to 0(e5)]: 

* +v .») Bp+ (e + vxB-^).^ = 0, («) 

in which the v x B term does not contribute since (v x B) ■ v = 0. The distribution function 
n p is decomposed as 

n p = + nf) + n^ 5) + ■■■ , (93) 

where 

f = Ofa - e p ) ~ OQi - ]p|) + - |p|), (94) 

which follows from the dispersion relation (43). Note that the second term in Eq. (94) is 
also 0(e5), but this is separated from n p ^ in our definition. In the calculations below, we 
have to add both contributions at the same order of 0(e5) [see Eq. (96) below]. 
The kinetic equation can be written down at each order as 



+ v -7r n p ) = E-v5(/i-|p|), (95) 



d d 
di + Y "dx 

1-1) {<" + 1> " IPl)) - v. £ (^) * - W) = 0. (96) 
The second equation is further simplified to 

^-5)"P--i(l?Kw>- (-) 

Using the method of characteristics, we can solve these equations; the operator in the 
left-hand sides of equations, v ■ d x , is the time derivative along the characteristic v = dx/dt. 
The solutions are given by 

POO 

= 5(fi- |p|) / dr e'^v ■ E(x - vr), (98) 
Jo 

poo 1 

n[f) = -6(fi- |p|) / dre^-vB^-DT), (99) 
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where 77 is a small positive parameter which ensures E(£ — > — 00, x) — > and B(t — > 
—00, x) — > 0. 

Now let us compute the current defined in Eq. (18). The current can be written down at 
the orders of 0(e) and 0(e5) respectively as 



Me) 



(x) 



d 3 p 

w. 

d 3 p 
(27p 



7; /v,( £ ) 



B i 



1 '. ; <;, '' H — |p|) — e 



ijk vi dn v ] 
2fi dx k 



(100) 
(101) 



where v 11 = (1, v). The zeroth component of the four-current of order 0(eS), i.e., the number 
density n {tS) (x), is found to vanish after the angular integration. 

First consider the four-current j^ e \ Substituting Eq. (98) into Eq. (100), the current 
reads 

d 3 p 



x) = 



(27T) 



gV^C/x — |p|) / dr e- 1,T v -E(x -vt). 



Using the useful formula 



J d 4 x e ik ' x 1^ 



dr e T,T f(x — vt) = 



if(k) 



(102) 



(103) 



v ■ k + IT) ' 

and the linear response theory (91), we obtain the retarded parity-even polarization tensor 
[40]: 



nf (fc) 



2tt 2 



5^5 



(104) 



47r v ■ k + ie 

which agrees with Eq. (84) derived from the perturbation theory. 

We then turn to the subleading three-current j l ^ 5 \ Substituting Eqs. (98) and (99) into 
(101), and using the formula (103), the current reads 



■i(e8) 



(k) 



d 3 p 



6(ii-\p\) 



ie klm v i v k ujk l A m ie ijk k j A k ie ikl v h k l [w(v • A) - (v • k)A ] 



2fx(v ■ k + ie) 



2/x 



2/j,(v ■ k + ie) 



(105) 



among which A term in the brackets vanishes after the angular integration. Using the linear 
response theory (91), we obtain the parity-odd polarization tensor 



TVi(k) 



An 2 



ie ijk k k + iu J 



dv (e jkl v l — e lkl v j )v k k l 



After the angular integration, this reduces to 



m(k) 



4tt 2 



-ie 



ijk j^h 



v ■ k + ie 



[l-uL(k)} : 



(106) 



(107) 



where L(k) is defined in Eq. (88). This is equivalent to Eq. (87) derived from the per- 
turbation theory, which confirms that the physics in the next-to-leading-order hard dense 
loop approximation can be described by the kinetic theory with Berry curvature corrections. 
Note that contributions of the inhomogeneous term in the current (21) and the magnetic 
moment in Eq. (43) are necessary for the matching of the correlation functions. 



19 



V. CONCLUSION 

In this paper, we have shown a way to bridge between quantum field theories and the 
kinetic theory with Berry curvature corrections that exhibits triangle anomalies and the 
chiral magnetic effect. The field theoretic procedure to derive such a kinetic theory devel- 
oped in this paper can, in principle, be generalized to higher order in gauge fields and/or 
derivatives. We have also computed the parity-odd correlation function using this kinetic 
theory, which is found to agree with the perturbative result beyond the leading-order hard 
dense loop approximation. 

It should be remarked that our derivation of the kinetic theory from underlying quan- 
tum field theories is limited to low temperature region T <C \x where the Fermi surface is 
well-defined and the high density effective theory is applicable; a generalization to higher 
temperature regime would be desirable. We also remark that our formulation based on the 
Berry curvature is not manifestly Lorentz covariant by construction. It might be possible 
to formulate the kinetic theory in a Lorentz covariant way similarly to the usual Vlasov 
equation. Without referring to the Berry curvatures, such a solution of the kinetic equation 
in the hydrodynamic regime is obtained in Ref. [11] which also reproduces the chiral vortical 
effect [20, 42-45]. 

The inclusion of collision terms in the kinetic theory (15) is straightforward, which gives 
the modified Boltzmann equation taking into account the anomalous effects. We hope that 
our work motivates numerical applications of the new kinetic equation (15) with or without 
collisions in various systems, such as dense quark matter, neutrino gas, early Universe at 
large lepton chemical potential, and doped Weyl semimetals. 
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